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Abstract. We describe a possible implementation of the nanomechanical quantum 
superposition generation and detection scheme described in the preceding, companion 
paper [Armour A D and Blencowe M P 2008 JVew. J. Phys. XX XXX]. The 
implementation is based on the circuit quantum electrodynamics (QED) set-up, with 
the addition of a mechanical degree of freedom formed out of a suspended, doubly- 
clamped segment of the superconducting loop of a dc SQUID located directly opposite 
the centre conductor of a coplanar waveguide (CPW). The relative merits of two 
SQUID based qubit realizations are addressed, in particular a capacitively coupled 
charge qubit and inductively coupled flux qubit. It is found that both realizations 
are equally promising, with comparable qubit-mechanical resonator mode as well as 
qubit-microwave resonator mode coupling strengths. 
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1. Introduction 

In reference [T], we described a scheme for generating and detecting superposition 
states of a nanomechanical resonator. The scheme employs a qubit with dispersive 
coupling to the mechanical resonator. By preparing the qubit in a superposition of 
energy eigenstates and the mechanical resonator in a displaced thermal state, the latter 
subsequently evolves into a superposition state, the existence of which can be inferred 
through appropriate qubit state control and measurement. 

In the present paper, we describe a possible implementation of this scheme (see 
figure CD) that is based on the demonstrated circuit quantum electrodynamics (QED) 
set-up [21 [3] . Control and readout of the qubit is achieved by pumping and probing the 
appropriate mode of a microwave coplanar waveguide (CPW) resonator that couples 
electromagnetically to the qubit. We shall consider two types of superconducting qubit: 
a capacitively coupled charge qubit, otherwise known as a Cooper pair box (CPB) [HE], 
and an inductively coupled flux qubit [6J. As we shall see, both types of qubit appear 
equally promising for implementation, having comparable strength couplings to the 
CPW resonator and to the nanomechanical resonator. One motivation for adopting 
the circuit QED set-up is the adequate demonstrated qubit coherence times for our 
purpose pp. 

We suppose that a segment of the qubit flux loop directly opposite the centre 
conductor of the CPW is freely suspended, forming a doubly-clamped beam mechanical 
resonator. Related schemes can be found in references [3 El El HDl E] • Furthermore, in 
addition to the usual GHz microwave qubit control/readout line that is capacitively 
coupled to the CPW, we require a separate, low frequency (i.e., tens of MHz and 
below) centre conductor bias line, in order to achieve the necessary strong couplings 
between the mechanical resonator and both qubit and CPW mode, as well as to drive 
the mechanical resonator on resonance, producing the displaced thermal state. In the 
case of the capacitively coupled charge qubit, the required centre conductor voltage 
bias can be introduced via a series inductor, otherwise known as a 'bias tee' that is 
microwave engineered to reflect the separate, much higher frequency GHz qubit control 
signal [12]. In the case of the inductively coupled flux qubit, the required current bias 
can be introduced via two superconducting low pass filters at each end of the centre 
conductor [13j . The high frequency thermal Johnson noise entering the bias line can be 
significantly suppressed with additional low-pass filtering at each temperature stage of 
the set-up, while at the same time allowing sufficiently large dc voltage or current biases 
for the required mechanical resonator-qubit coupling strengths [14] . On the other hand, 
the high frequency thermal noise entering the microwave control line can be reduced 
by attenuating the signal [H]; only low drive power signals are required to control the 
qubit 0. 

Sections 2 and 3 of this paper derive the respective charge and flux qubit 
Hamiltonians with coupled, single microwave and mechanical modes. Both qubit 
Hamiltonians have a common origin in the CPW-coupled-dc SQUID circuit of figure [TJ 
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Figure 1. Model circuit of the closed CPW-dc SQUID system. The CPW comprises 
a uniform width centre conductor of length I and parallel ground planes. The SQUID 
is located at the midway point (z — 0) in the gap between the centre conductor 
and one of the ground planes. The mechanical resonator is formed out of a freely- 
suspended segment of the SQUID loop located directly opposite the centre conductor. 
An electromagnetic mode of the centre conductor couples both capacitively (via mutual 
capacitance C m ) and inductively (via mutual inductance M m ) to the SQUID loop 
circulating currents. Not shown are the microwave qubit control/readout line and low 
frequency centre conductor bias line. 

Section 2 discusses the classical dynamics, while section 3 derives the quantum 
Hamiltonians of the charge and flux qubit systems, equations ( 1351) and ( j4"6j) respectively. 
While there is nothing particularly new in terms of theory in these two sections, they 
serve the purpose of mapping the scheme of reference p] onto the circuit QED set- 
up with a unified treatment of the charge and flux qubit implementations, facilitating a 
direct comparison between them. Readers not interested in the derivations may proceed 
directly to section 4, where the relative merits of the qubit Hamiltonians for actual device 
implementation are addressed. Section 5 gives a brief conclusion. Appendix A gives a 
derivation of the single mode approximation to the CPW superconducting phase wave 
equation, while appendix B derives the qubit Langevin equations. 

2. CPW-dc SQUID classical equations of motion 

In this section, we analyze the classical equations of motion for the phase field 
coordinate <fr(z,t) of a length I coplanar waveguide (CPW) with longitudinal coordinate 
—1/2 < z < 1/2 and gauge invariant phases 4>i an d 4>2 across each of the two Josephson 
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junctions (JJ) of the dc SQUID. These equations are the common starting point for the 
subsequent derivations of both the microwave resonator mode-flux qubit and microwave 
resonator mode-charge qubit Hamiltonians. Related analyses can be found in the review 
references [151 EE El EH] • The CPW and SQUID are coupled via a mechanical, position- 
dependent mutual inductance M m and capacitance C m at z — 0, arising from the 
adjacent, parallel centre conductor and suspended SQUID loop segments. The position 
dependence and hence mechanical degree of freedom will be suppressed for the time 
being, to be introduced in section [3] once the qubit Hamiltonians have been obtained. 
Figured] details the closed CPW-SQUID circuit. Not shown are the capacitively-coupled 
microwave signal pump-probe line, as well as the low frequency bias line. Again, we 
will ignore these external lines, focusing first on the closed CPW-SQUID dynamics. 
The external bias dependences will be introduced later below once the relevant classical 
CPW-coupled SQUID Lagrangians have been derived [see equations ( 1T6]) and (1271)]. 

For — 1/2 < z < and < z < 1/2, the CPW phase field coordinate obeys the wave 
equation 

dt 2 V dz 2 ' 1 ' 

where v = (L W C W ) -1 / 2 is the phase velocity and L w and C w are the inductance and 
capacitance per unit length of the CPW, respectively. At the CPW boundaries z = ±1/2, 
the current /„, = — „ vanishes, since as mentioned above we are in the first instance 
considering the CPW-SQUID system to be closed. We assume small CPW-SQUID 
couplings, i.e., C m /(C w l) <C 1 and M m /(L w l) <ti 1. 

The equations of motion for the coupled SQUID can be conveniently derived by 
applying Kirchhoff's laws for the voltages and currents at z = ± and in the SQUID 
loop. Referring to figure [□, current conservation gives 

J--/+ = / 1+ / 2 . (2) 

The voltage relations are 

v _ = M m dl x Q m M m dl~ Ldh $ #i 

w ' 2 dt C m 2 dt 2 dt 2?r dt U 

and 

. = _Mrndh Qm MradJ^ Ldh £o d<fo 

2 dt C m 2 dt 2 dt 2tt dt ' U 
where V w = f^f^, Qm is the charge on the capacitance C m and L is the SQUID loop 
self inductance. Current conservation at the JJ nodes gives 

h = i c sin 0i + — — 0i 

c % 

I 2 = I c sin 02 + —^f- '^ ( 5 ) 

where I c is the JJ critical current and Cj the JJ capacitance (assumed identical for each 
junction). The J J phases are related to the net flux $ threading the SQUID loop as 

2 - 0i = 2vm + 27T-J-, (6) 

$0 
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where the net flux is 

$ = $ext 



M„ 



(7) 



with $ext the externally applied flux. It is more convenient to work with the SQUID 
phase variables 7± = (0i ± 02) /2, in which case equation (j3J) becomes 

h + h I--I+ 



and 



uij 7+ + sin 7 + cos 7_ 



to j 7_ + cos 7 + sin 7. 



2/, 2/, 



= M m I- + I+ 
L 2I C 

and ([ZD, = 



2 



7- + ir n + 



(8) 



(9) 



where we have used equations flBj) and (JTj), u;j = ^/2ttI c / (Cj$ ) is the J J plasma 
frequency and = 2tiLI c /$q. Subtracting and adding the voltage relations and 
(jlj) in order to express them in terms of variables j±, we obtain 



V~ - U H 



M 2 

m 

2L 



C + It + 



M m $ d 



7- + 7T n + 



(10) 



where we have again used equations 

1 



V~ + 

r 7/) 1 r 11} 



(L — 2M m ) I 



-7+, 



L Ti dt 
and ((Tj), and 

where we have taken time derivatives to replace the charges with currents and we have 
also used the current conservation relation (j2J). Equations (|HJ), ©, ({10]) . (TlTI) along with 
the wave equation (CD) completely specify the classical dynamics of the coupled CPW-dc 
SQUID system. 

We shall restrict ourselves to CPW mode solutions that reflect the symmetry of the 
circuit. In particular, we consider two types of solution: 1) Voltage antinode and current 
node at z = 0: V~ = V^~, I~ = —I^\ 2) Voltage node and current antinode at z = 0: 
V~ = —Vj~, I~ — I+. In the former case, the CPW-SQUID coupling is predominantly 
capacitive and results in the microwave resonator mode-charge qubit Hamiltonian, while 
in the latter case the coupling is inductive and results in the microwave resonator mode- 
flux qubit Hamiltonian. The following two subsections deal with each case in turn. 



2.1. Voltage antinode equations: capacitive coupling 

For a voltage antinode-current node solution, we have from equation ( flQl) that 7_ must 
be a constant. From equation dHJ), we then require that 7_ = qn and $ ext = — (g + n)$ , 
where q is an arbitrary integer. Thus, we see that in order to have a non trivial solution 
for the remaining 7+ variable, the external flux is constrained to be an integer multiple of 
the flux quantum. However, if (3l <C 1, i.e., the SQUID loop self inductance L is small, 
then to leading, zeroth order in (3l a voltage antinode solution is allowed provided 
7_ = — 71 (n + $ ext /$o) an d the external flux need not be constrained. Assuming 
therefore a small loop area SQUID with (3l <C 1 in the CPW resonator voltage antinode 
case and neglecting the inductance dependent terms in equation pip , we obtain 

V- = ^-1- + %+■ (12) 
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for 7+ thus becomes 
'-Ct,1+ + 2I C cos (7r$ ext / $ ) sin 7+ 



C m V w , 



(13) 



Using equation (TT21) . equation 
2tt 

where C*s = 2Cj + C m and we have set n — 0, since observable quantities do not depend 
on n. Equations (|T2|) and (|T3|) along with the wave equation ([1]) completely specify the 
voltage antinode restricted classical dynamics. 

The qubit control and readout protocols [1] require driving the CPW close to one of 
its resonant modes. In the voltage antinode case, this is usually the second fundamental 
mode [3]. The subsequent analysis of the CPW-SQUID dynamics is greatly simplified if 
we can first effectively replace the CPW with a single harmonic oscillator at the relevant 



microwave mode frequency. In |Appendix A[ we derive these approximate, effective single 
oscillator mode equations and also give the conditions under which these equations 
are valid. In the voltage antinode second fundamental mode case, we obtain for the 
approximate coupled single oscillator-SQUID equations of motion: 



Cs^7+ + 2/ c cos(7r$ cxt /$ ) sin 7+ = C r 

Z7T 



$ 

'27 



and 



CJ $ 



o 



2 2tt 



+ 



2ttv 



n $o.. 

C m— 7+) 
Z7T 



(14) 



(15) 



where <p(t) is the oscillator coordinate giving the CPW phase field amplitude at the 
SQUID location z = and where we have neglected an order C m /(C w l) <C 1 shift in 
the CPW second fundamental mode frequency uj w = 2ttv/1. 

Equations of motion (TBI) and ffTBI) follow from the Lagrangian 

2 



£(7h 



■ 7+, 



v r° 



7T 



$0 

2tt 



cos(7r$ext/$o) cos 7+ 



2^^ (j) 

T 



2tt ( 



+ V(t) 



(16) 



where we have included a CPW voltage bias V(t), assumed to be slowly varying in 
time as compared with the microwave mode frequency uo w . As we shall see in section HI 
this additional voltage bias is necessary in order to tune the qubit 's operating point, 
couple strongly the qubit with the mechanical oscillator, and also drive the mechanical 
oscillator on resonance pp. 

Using the Lagrangian ( fTBl to construct the associated Hamiltonian, we have for the 
generalized momenta: 



and 



V<i> 



P74 



dL _ 
dL 



C v , I 



£0 
2tt 



£0 

2tt 



Cr, ^ 



$0 
2tt 



7h 



r 



£0 

2tt 



7+ 



Lix 



(17) 



(18) 
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Expressing the velocities in terms of the momenta and coordinate variables and 
substituting into the definition for the Hamiltonian, H = p^tfi + p 7+7+ — L, we obtain 
after some algebra 

H{7+,<t>,P 7+ ,P<t,) = Cj, 1 f^J ^ - : ^cos(7r$ cxt /$ )cos7 + 



CJ \%J 2 2 V 27r / V I 



Cm. ( 27T 



+ 2 cjci{rJ < 19) 

where we have have dropped overall constant (i.e., variable independent) terms. 

From the definition (|18p for the generalized momentum p 1+ and recalling that 
$o = hj (2e), we see that 

P 1+ = -H[N - N m (t)], (20) 

where N is the number of excess Cooper pairs on the island represented by the 
directly connected nodes between the C m and Cj capacitances (see figure [T]) and 
N m (t) = —C m V(t)/(2e) is the polarization charge induced by the slowly varying voltage 
bias on the C m electrode, expressed in units of Cooper pair number. Expression (T20l) 
then suggests an alternative formulation of the Hamiltonian in terms of the Cooper pair 
number N: 

if ( 7+ , 0, N, P(f> ) = 2$-[N- N m (t)} 2 - — cos(7T$ cxt /$ ) cos 7+ 

2C E 7T 

2 {2n\ 2 pl 1 , /$ n \ 2 /27tu n 2 a2 



CJ \$ J 2 2 V2vr / 
-4e^[^V-iV m (t)]- P , (21) 

and we recognize in the top line of equation (T2"Tj) the usual Cooper pair box (CPB) 
Hamiltonian II!). 



2.2. Voltage node equations: inductive coupling 

For a voltage node-current antinode solution, we see from equations (TTTT) and (jHJ) that 
7+ = g7r, where g is an arbitrary integer. Thus, equation (Q for 7 _ becomes 

-2- , • 2 / 7T$ cxt \ M m /~ 



^_ + sin7 _ + _^_ + ^. Tt , (22) 

where the mutual capacitance C m drops out and we have set n = q = 0. Equation (fTOl) 
also simplifies to 

K, = " 2T/« + T ^T-, (23) 

where we have assumed that $ ext ~ 0, i.e., changes slowly compared with the dynamical 
timescales of the SQUID and CPW. Equations (1221) and (1231) along with wave equation 
(TTJ completely specify the voltage node restricted classical dynamics. 
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The approximate, first fundamental (i.e., lowest frequency) microwave mode- 
coupled SQUID equations that follow from equations (1221) and (1231) are 
(see l Appendix A ) : 

2 



sin 7_ + 



0L 



7- + 



7T$ r 



o 



L 2ttL, 



and 



4 2tt 



tan(fc // 2 ) 



/7TU V 

+ (t) 



o 



L 2ttL, 



■&o7- 



■fcotan(fc o J/ 2 )0( 24 ) 



(25) 



where from equation ( 1A.4I) in Appendix A the term $o/(2vrL^)fco tan(fcoZ/2)</> appearing 
in the 7_ equation of motion (124"j) is just the CPW current at the SQUID location z — 0. 
The wavenumber fco is the fundamental solution to the equation 



tan 



2 



LyjlL 



2 — 12/ M 2 (26) 

Neglecting order M m j (L w l) <C 1 corrections, we have approximately for the fundamental 
wavenumber and mode frequency: k = ir/l and u w = ttv/1, respectively. 
Equations of motion (1241) and (1251) follow from the Lagrangian 



/$ y. 2 j c $ x ( 



2tt 



H 7- 

L 2n 



2n 

tan 2 (fc //2) 
$ 



7T 



[$ ext - M m /(t)] 



/7TV\ 2 (fr 



2ttL„ 



fc tan(A; Z/2)</> 



(27) 



where we have included a CPW current bias /(£), assumed to be slowly varying in time 
as compared with the fundamental microwave mode frequency uj w . As we shall see in 
section HJ this additional current bias is necessary in order couple strongly the qubit 
with the mechanical oscillator and also drive the mechanical oscillator on resonance pQ. 
The Hamiltonian associated with Lagrangian ([271) is 



'2tt\ 
$ J ^ 



+ 



4 



C w l 
M m $ 



2tt 

2 (A,//2) ( — 
$ 



-i 



tan-^oV 2 ) 
$ 



cos7_ - (3 L 
2 + 4 J 



71 



+ [$ext - M m I(t)\ 
<P 



$0 

2tt 



tan 2 (fc //2) 



TTV\ 2 (f) 2 

t) t 



L 2tt 



7- 



2ttL, 



/cotan(fc //2)</) 



(28) 



3. Quantum equations of motion 



The goal of this section is to derive the quantum Langevin equations for the coupled 
CPW microwave oscillator-charge(fTux) qubit-mechanical oscillator system, taking into 
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account the effects of damping and noise on the three subsystems. We shall begin 
by deriving the closed system quantum Hamiltonians, treating the (voltage antinode) 
capacitively-coupled charge qubit and (voltage node) inductively-coupled flux qubit 
systems separately as we have done in the previous section. 




3.1. Charge qubit 

It is most convenient to work in terms of the CPW microwave oscillator raising and 
lowering operators: 

^ = /o i ( m ^ u ^ =F Wt) > ( 29 ) 
V lm w u w n v / 

where, from equation (12ip . m w = C w I/2(^q/2tc) 2 and u w = 2nv/l. The Hamiltonian 
f l2Tj) then becomes 

H = { ^r\ N ' N m{t)] 2 - — cos(7r$ cxt /$ ) cos 7+ + hw w a+a w 

71 

(2e) 2 

±-L[N-N m (t)}iK-a w ), (30) 

where for notational convenience we have dropped the hats and also the '-' superscript on 
the lowering operator. Examining the interaction Hamiltonian, we see that the coupling 
constant in units of Huj w depends essentially on the ratio of the CPW single Cooper pair 
charging energy to the CPW mode single photon energy [3]. 

Note from equation (1201) that the standard Poisson bracket relation { 7+ ,P 7+ } = 1 
gives { 7+ ,iV} = —fr 1 for the variables 7+ and N. Thus, when we quantize using the 
correspondence principle between Poisson brackets and commutators, i.e., {•,•} = c — > 
— iftr l \-, ■] = c (for some constant c), we therefore have that [7+, N] = —i. 

We now truncate the SQUID Hilbert space down to the lowest energy, two- 
dimensional Hilbert space, obtaining the charge qubit Hamiltonian. Working in 
the representation for which the number operator N is diagonal, the commutation 
relation [ 7+ ,iV] = —i gives 7+ <-> —id/dN. For the number operator, we have 
N = N\N)(N\ + (N + 1)\N + 1)(N + 1\ = (N + |)(|iV)(JV| + \N + 1)(N + 1|) - 
l(\N)(N\ - \N + 1)(N + 1|) = (N + \)I - \o z , where we have truncated to the 
two-dimensional Hilbert space spanned by eigenkets \N) = ( ), \N + 1) = (°) 
of the number operator; in order to be within this subspace, we require V(t) to 
be tuned such that N < N m (t) < N + 1, where recall N m (t) = -C m V(t)/(2e). 
The truncated charging energy term can similarly be expressed as [N — N m (t)} 2 = 
(N - N m ) 2 \N)(N\ + (N + 1 — N m ) 2 \N + 1)(N+ 1|) = [5N{t)] 2 + \ + 5N(t)a z , where 
SN(t) = N m {t) - (JV+ I). Writing cos( 7+ ) = (e^+ + e~^)/2 = {e d ' dN + e~ d / dN )/2, we 
have cos( 7+ ) = (\N)(N + 1| + |JV + l)(JV|)/2 = o~ x /2. Substituting the above truncated 
approximations into equation fl3"Ul) . we obtain for the truncated Hamiltonian: 

H = E Cs 5N(t)a z - ^Ej($ cxt )o- x + kuj w a+a w + hg c {a w + a+)a z 

+ 2hg c (a w + a+)SN{t) + E Cs [5N(t)] 2 , (31) 
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where 

C m l (2ey/(CJ) 

9c = Uw 2^\j ET W (32) 

is the microwave mode-qubit coupling [5], and where E Ct . = (2e) 2 /(2Cs) is the 
single Cooper pair charging energy for the total capacitance = 2Cj + C m and 
Ej($ext) = 2Ej cos(7r$ ex t/$o) is the effective flux-dependent Josephson energy, with 
Ej = J c $ /(27t) the Josephson energy of a single J J. Note that we have also made the 
canonical replacement 

We now incorporate the mechanical degree of freedom, supposing that a segment of 
the SQUID loop directly opposite the centre conductor forms a mechanically compliant, 
doubly-clamped beam. For a small in plane displacement with centre-of-mass amplitude 
x, the position dependent mutual capacitance is C m (d + x) ~ C m (d)(l — r/cx/d), where 
d is the static equilibrium gap between the centre conductor and opposite facing SQUID 
loop segment and 7]c is a geometrical factor of order one accounting for the non-uniform 
beam displacement. Substituting this C m (d+x) into the Hamiltonian (I3T]) and assuming 
that C m <C C*e, we see that to leading order the coupling to the mechanical resonator 
arises through a modulation of the V(t) induced polarization N m (t), as well as through a 
modulation of the microwave mode-qubit coupling g c . Express the applied voltage bias 
as V{t) = Vdc + V ac (t), where Vd c is a constant dc voltage and V ac (t) is slowly varying 
as compared with the microwave mode frequency uo w . We obtain for the mechanical 
resonator contribution to the Hamiltonian: 

H m = huimO^am + X c (a m + a+)cr 2 - hx c gc(a w + a+)(a m + a^)a z 

, C rn V ac (t) . i , /oo\ 

-Ac (a m + <J, (33) 

e 

where xq = Vc x z P /d with x zp the zero-point uncertainty and we assume that it 
is the frequency u m fundamental flexural mode of the mechanical resonator that 
predominantly couples to the qubit with strength 

Ac = Vc^^eV dc . (34) 



Hamiltonian (1331) neglects direct microwave-mechanical oscillator interaction terms, 
assumed to have a small effect on the coupled dynamics since the two modes are 
significantly off- resonance: u> m u w . We also assume that \Vd c \ 3> e/C m , necessary in 
order to achieve strong mechanical oscillator-qubit coupling, and that the mechanical 
oscillator driving voltage |V^ C (£)| e/C m . Later, we shall see that the typical large 
mechanical quality factors ensure that such ac driving voltages are adequate for exciting 
the mechanical oscillator to sufficiently large amplitude as required by the qubit control 
protocol p. 

The control protocol also requires operating at the charge degeneracy point SN = 
—C m Vdc/{2e) — N = 0, where the CPB qubit coherence time is a maximum [20]. Thus, 
the full microwave-qubit-mechanical oscillator Hamiltonian becomes 

H = \^w a a z + hu w a+a w + hw m a+ l a m + hg c {a w + a+)a x 
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+ Ac(a m + a+)a x - hx c g c (a w + a+)(a m + a^)a x 
, C m V ac (t) 

-Ac (a m + a„), (35) 

where we have changed the Pauli matrix basis: a x — > — cr 2 , cr 2 — > — a x , and cj a = 
Ej($ ext )/H is the qubit transition frequency. From (1351) . we see that the mechanical 
oscillator couples to the qubit and microwave oscillator through two interaction terms of 
the form {a m + a^)a x and (a w + a+)(a m + a+jaj., respectively. Given the wide separation 
of timescales, u m Cw ffl , an adiabatic approximation can be performed to show that the 
((tmlo^ffj interaction term gives rise to an approximate, dispersive interaction a+a m <7 2 
between the mechanical oscillator and qubit [7j. Thus, a dominant effect of the qubit 
on the mechanical oscillator is to shift its frequency up or down depending on the state 
of the qubit. By preparing the qubit in a superposition state, the mechanical oscillator 
is in turn driven into a superposition state via the dispersive interaction, as analyzed 
in reference pQ. Provided the frequency difference \u w — u a \ is much larger than the 
coupling ga, the microwave-qubit interaction (a w + a+)a x can similarly be replaced by 
an approximate dispersive interaction a+a w a z j3]. Supposing that the qubit remains in 
its ground state, the microwave mode then exerts a ponderomotive-type force on the 
mechanical oscillator described by the interaction a+a^am + a^). Thus, passive cooling 
or alternatively amplification of the mechanical motion can in principle be implemented 
with appropriate red or blue detunings of the microwave drive frequency with respect 
to the microwave mode frequency u w [ID] . We discuss these possibilities in section HI 

3.2. Flux qubit 

In this section, we obtain the truncated flux qubit Hamiltonian. While the analysis 
is well known and involves basic Schrodinger wave mechanics, it does serve a purpose 
in bringing to the fore the differences in the flux and charge qubit parameter regimes. 
Readers not interested in the details of the derivation may skip directly to the resulting 
qubit Hamiltonian (|46|) . 

From equation (|28|) . the microwave mode effective mass is 

m w = -CJ (—} tan 2 (A;o//2) (36) 



4 \2ix 



and the microwave mode fundamental frequency is u m = nv/l. In terms of the raising 
and lowering operators, the Hamiltonian (l28j) is then 



H = Cf 



2vr \ 2 V%_ i c $c 



$n 7 2 2tt 



. 2' 

cos 7 _ - Pi 1 ( 7 _ + [$ ext - M m I{t)} 
<P 



+ huj w aZa w —d2hhJ w — — — 7_(a^ + a tu ). (37) 

Examining the interaction Hamiltonian, we see that the coupling constant in units of 
h\o w depends essentially on the ratio of the CPW single flux quantum addition energy 
to the CPW mode single photon energy. 
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The eigenstates and eigenvalues of the SQUID part of the Hamiltonian (137)) can be 
obtained, e.g., by working in the Schrodinger wavefunction form. In the representation 
7^7 and p 7 <-> —ihd/d'y (where we have dropped the ' — ' subscript on 7 for notational 
convenience), the SQUID Schrodinger equation becomes 



di 2 



+ Ej 



Pi 



TV 



7 H [$ cxt 

' $0 cxt 



M m I, 



cos 7 



¥( 7 ) 



= ^(7), (38) 

where E Cj = (2e) 2 /(2Cj) is the Cooper pair charging energy of the JJ capacitor, 
Ej = J c $ /(2e) is the Josephson energy, and we have suppressed for the present the 
slowly time- varying contribution I ac (t) to the current bias 1(f) = I^c + I a c(f)- Provided 
(3 L = 2tcLI c /& > 2 and $ ext -M m J dc = (2n+l)$ 0) n = Q, ±1, ±2, then the potential 
energy function admits a symmetric double well centered at 7 = — (2n + 1)tt. Figure [2] 
shows an example plot of the potential V(7) in region of its minimum, indicating a 
double-well. Also shown is the resulting asymmetric double well when the external flux 
$ ex t is detuned slightly from the above symmetric well condition value. In order that the 




Figure 2. Double well potential for (3l 
$cxt - M m I dc = -0.97$ (dashed line). 



3 and <i> c 



M m I d 



-5>o (solid line), 



lowest energy eigenvalues lie 'within' the double well and hence be strongly anharmonic, 
we require that the harmonic ground state energy E neglecting the cosine term in the 
potential V(7) be smaller than the central barrier maximum (= Ej for the degenerate, 
symmetric well — see figure |2J). Thus, we require that E = a/ 'E Cj Ej / '/3 L < Ej, or 
Eqj/Ej <C However, Ec, cannot be too small, for otherwise the energy spacing 
between the ground and first excited states will correspondingly be too small, making 
the resulting qubit susceptible to damping and decoherence by thermal and other low 
energy noise sources. 

We have seen that, in contrast to the charge qubit (see section 2.1), the flux qubit 
requires a nonnegligible @l, i>e, the flux qubit requires a much larger SQUID loop self- 
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inductance L than the CPB. Adequate effective self-inductances can be achieved by 
using small loop, multiple Josephson Junction SQUID devices [6], where the additional 
JJ's compensate for the small loop area, hence avoiding potential increases in flux noise 
that come with using larger loop areas. However, in the present work, we consider only 
a double JJ, dc SQUID for simplicity; the relevant multiple JJ SQUID devices do not 
involve any conceptually new features and can be analyzed along similar lines. 

It is quite informative to solve the Schrodinger equation (1381) approximately in terms 
of the two harmonic ground basis functions for the quadratic well expansions about the 
two minima. We have: 



^(7) = V( % ) + Ej\(5l l + cos(7,)/2]( 7 - 7 *) 2 , 
where the double well minima, ji, i = 1,2, are solutions to 



20i 



71 



7i + — [$ext - M m i d< :\ 

<P 



+ sin 7i = 0. 



(39) 



(40) 



The normalized ground eigenfunctions of the harmonic approximation Hamiltonians 
with V( 7 ) replaced by V(7i) are 



Ml) 



E, 



TT 2 E, 



Cj 



(^ 1 + cos( 7 ,.)/2) 



1/8 



x exp 




[^ 1 + cos(7,)/2]( 7 -7,) 2 



(41) 
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B 








The lowest two eigenfunctions of the double well Hamiltonian H are then approximately 
expressed as linear combinations of l^j): |^) = ^l^i) + i? | ^2) , where the constants A, 
B are solutions to 



(42) 



and the solutions E give the approximate energy eigenvalues of H. Note that the 
two basis functions are not orthogonal; the solutions are a reasonable approximation 
provided |(^i|^ 2 )| < 1- 

Let us now consider a concrete example. Taking L = 500 pH and I c = 2.0 fiA, we 
have Pl = 3.03 and Ej = 4.1 meV. Taking Cj — 1 fF, we have E Cj = 0.3 meV, so that 
E Cj /(Ejf3 L ) = 0.03 < I. Considering a symmetric well with $ ext -M m J dc = (2n+l)$o, 
we find that A = ±5, and the approximate energy eigenvalues are: £L = 3.95 meV and 
E + = 4.01 meV with difference AE = E + - E_ = 55 /xeV (and AE/h = 13 GHz). The 
approximate normalized eigenstates are \^f±) = 0.71(|\E f 1) ± 1^2)); i>e., (anti) symmetric 
combinations of the basis states, with the symmetric combination associated with the 
lower, ground eigenvalue E_. For this example, we have (^1^2) — 0.01. On the 
other hand, for a slight external flux bias detuning $ ext — M m I dc = (2n — 0.98)$ ) 
making the double well potential asymmetrical, we have E_ = 3.75 meV and E + = 
4.20 meV with difference AE = E + — E- = 450 /xeV, and associated eigenstates 
= 0.998|^i) - 0.058|^ 2 ) and |*+) = 0.07|^i) + 0.998|^ 2 )- Thus, the lowest 



Probing the quantum coherence of a nanomechanical resonator II: Implementation 14 



two eigenstates correspond to the phase being approximately localized to one or the 
other of the two wells. 

Proceeding formally with the truncation of the Hamiltonian ( 1371) to the subspace 
spanned by exact eigenstates \^±) and incorporating the mechanical degree of freedom 
through a position-dependent mutual inductance M m (d + x) ~ M m {d){l — T]Lx/d) (with 
t]l a geometrical factor of order one), we have: 

H = -AEa z + fruj w ala w + fko m a^a m - hg L (a w + 0^)7 

+ r) L —f———I dc (a m + a£) 7 + -H$cxt - M m I dc \ 
a L Zn \ <Po / 

r M 

+ Vl—, — j^hg L (a w + aZ){a m + 0^)7 

+ VL^f^^Iac(t)(a m + o+) (7 + - 2M m I dc \) , (43) 

a L Zn \ <Po / 

where the microwave mode-qubit coupling is 

M m / 2($ /27t)V(LJ)" 

^ = W "-TV ^ {) 

and the truncated phase operator 7 is 

7 = l+-<*x + ^(7++ - 7— + ^(7++ + 7— (45) 

with 7 ++ = (\I> + |7|\I/ + ), etc. Note that AE and the truncated matrix elements 7^ are 
$ ex t and /dc-dependent. For the example parameters considered above with symmetric 

double well potential [$ ex t — M m I dc = (2n + l)$o], we have: 7++ = 7 = — (2n + l)n, 

and 7 + _ ~ —n/2, so that the o z operator does not arise in the couplings. On the other 
hand, for the above asymmetric well example [$ ex t — M m I dc = (2n — 0.98)$o]> we have 

7 ++ = 4.6 — 2nn, r y = 1.6 — 2nn, and 7 + _ = —0.2, so that the a z operator term is 

now present with much larger coupling as compared with the a x operator term. 

We assume that the flux qubit's coherence time is optimized for the symmetric well 
case, i.e., the bias conditions are: $ e xt — M m I dc = (2n + l)$o- The Hamiltonian ( 1431 ) 
then simplifies: 

H = ^hw a a z + hw w a+a w + huj m a^a m - hg L ^ + _(a w + a+)a x 
+ A L 7+_(a m + a+)a x + hx L gLl+-{a w + <4>)( a m + a^)a x 

~ Xl (eh I \ ( a ^ + Q m) ; (46) 

where u a = AE/h, xl = f]LX zp /d, and the qubit-mechanical oscillator coupling is 

x zp M m $ I dc 

XL = r]L lT — ^2n- (47) 
In the Hamiltonian simplification (|46p . we neglect direct microwave oscillator-mechanical 
oscillator interaction terms and also assume small mechanical oscillator drive currents 
|i«c(t)| < $o/M m . 
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Note that the microwave oscillator-flux qubit-mechanical oscillator Hamilto- 
nian f|46|) is formally identical to the microwave oscillator-charge qubit-mechanical os- 
cillator Hamiltonian ( 1351) . the differences occuring only in their respective coupling 
strengths. This close identification will enable a straightforward comparison of the 
relative merits of the two Hamiltonians in section HI 



3. 3. Open system quantum equations of motion 

We now include model environments for the CPW-qubit-mechanical resonator system. 
For the CPW, we suppose that the dominant direct dissipation mechanism is due to 
coupling to the microwave pump-probe line used to drive and read out the state of 
the CPW. Assuming weak mechanical oscillator-bath and CPW mo de-pump /probe line 
couplings, with the baths modeled as a continuous spectrum of free harmonic oscillators, 
we obtain the following Langevin equations for the mechanical and CPW mode operators 
a m and a w , respectively: 

= ~[a m , H] - lm a m - i^2^e^a™ (48) 

and 



- '/,„, H] - lw a w - iy/Tn^e^cfc, (49) 



where H is the charge qubit Hamiltonian 035)) or flux qubit Hamiltonian (1461) . The 7j's 
are the various mode amplitude damping rates and the 'in' and 'out' bath operators are 
related as follows |23j : 



aT\t) - a?(t) = -i^2j i e- l<l>i a i (t). (50) 

This identity is used to determine the measured quantities' expectation values involving 
a° ut (t) in terms of the solutions to the mode operators cij(i) and the prescribed initial 

4°(t)' S . 

The qubits themselves also directly couple to environments other than the damped 
microwave and mechanical oscillators. For completeness, Langevin equations describing 



the dissipative qubit dynamics are given in Appendix B although they are not used in 
the present paper. 



4. Comparison of the flux and charge qubit implementations 

In the previous section, we derived Hamiltonians describing a charge qubit ( 1351) . or 
alternatively flux qubit (146]) . coupled to a microwave oscillator and mechanical oscillator. 
For convenience, we reproduce the two Hamiltonian expressions here: 

-^charge = ^hw a a z + hw w a^,a w + hiOraa^dm + hg c {a w + a^)cr x 
+ X c (a m + ati) a x - hx c gc(a w + a+)(a m + a^)a x 
- Ac [a m + a£) 
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and 

-ffflux = —hOaPz + hu w a+a w + tvjj m a+a m - hg L ^ + ^{a w + a+)a x 

+ A L 7 + _(a m + a^)a x + hx L g L -y + ^(a w + a+)(a m + a+)<r x 
M m I ac {t) t 

~ L ($oA) 

An important question concerns the relative merits of these two implementations. In this 
section, we address this by comparing realisable microwave-charge qubit and microwave- 
flux qubit coupling strengths go (1321) and g^ (JED, qubit-mechanical oscillator coupling 
strengths Ac (1341 and A^ (H7|) . and realisable mechanical displacement amplitudes due 
to achievable drive strengths V ac (t) and I ac {t)- 

Refering to the CPW circuit in figure [U we consider a geometry modeled on 
the circuit QED device of reference [21] dictated by the requirement that the CPW 



impedance Z = ^JL W /C W matches as closely as possible the typical 50fi impedances 
of coaxial feedlines. Furthermore, the resonant frequences of the first (voltage node at 
midpoint) and second (voltage antinode at midpoint) modes are required to be in the 
5 — 10 GHz range in order be in the quantum regime for the CPW at low dilution fridge 
temperatures, as well as to match the level separations of the flux or charge qubit. We 
consider the centre conductor of the CPW to be uniformly 10 /iiii wide and separated 
from the lateral ground planes by a 5 /zm gap. The doubly-clamped, suspended beam 
segment of the dc SQUID directly opposite the centre conductor is assumed to be 200 nm 
wide, with a gap of 100 nm between the two. All the described elements of the CPW 
and SQUID are 200 nm thick. The substrate is assumed to be silicon, with an etched 
out trench beneath the mechanical resonator that extends up to the adjacent centre 
conductor edge. 

The mutual capacitance per unit length between the centre conductor and lateral 
ground planes, C w , is obtained as follows. Apply two linearly independent static 
voltage configurations to the centre and lateral planes and for each configuration 
evaluate (numerically) the total stored electric field energy per unit length: U = 
\Ym j=iViC wi jVj. Next, invert to obtain the capacitance per unit length matrix C wi j. 
Finally, obtain the mutual capacitance per unit length in terms of the capacitance matrix 
using 

On the other hand, to obtain the self inductance per unit length of the centre inductor, 
L w , apply some non-zero steady state current / to the centre conductor, with zero 
applied current through the lateral conductors and then evaluate (numerically) the 
total stored magnetic field energy per unit length U = \L W I 2 . The self inductance 
is then simply L w = 2U/I 2 . Carrying out these procedures, we obtain C w = 
2.01 pF/cm and L w = 6.35 nH/cm. These values give an impedance Z = 56f2 and 
v = (LyjCw)' 1 ! 2 = 8.85 x 10 7 m/s for the phase velocity. Therefore, the first mode 
frequency is f\ = v/(2l) = 4.43 GHz/cm and second mode frequency is twice this: 
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/ 2 = v/l = 8.85 GHz/cm. Thus, the CPW should be about 1 cm long in order to have 
the first two microwave modes in the desired frequency range. In order to evaluate the 
mutual capacitance C m and mutual inductance M m between the centre conductor and 
suspended SQUID beam, we assume for simplicity an infinitely long beam and obtain: 
C m / 1 « 50 aF//zm and M m /l « 1 pH///m. 

For the microwave mode-qubit coupling strengths, we also require the SQUID qubit 
JJ capacitance values (for the antinode case) and self inductance values (for the node 
case). As example values, we refer to the CPB qubits of the Yale group [2] and the 
flux qubits of the Delft group [6] . We shall consider only order of magnitude estimates. 
Beginning with the inductively coupled flux qubit, we first note that the requirement 
Pl ~ 1 gives L ~ glr - Thus, in order that the SQUID inductance L and hence loop 
area not have to be too large, we require JJ's with large critical current I c values. The 
Delft group JJ's have I c ~ 1 /xA, giving L ~ 100 pH. Considering a CPW length 
of 1 cm gives a fundamental (node) mode frequency v/(2l) ~ 5 GHz, comparable to 
the Rabi oscillation frequency at the symmetry point of the flux qubit in reference [B|. 
A few micron long mechanical resonator gives M m ~ 1 pH, so that we have for the 
inductance ratio M m /L ~ 0.01. For the capacitive coupling in the case of the CPB 
qubit, a CPW length of 2 cm gives a second (antinode) mode frequency v/l ~ 5 GHz, 
comparable to the single electron charging energy (in units h~ r ) of the original CPB 
qubit in reference [2] . Therefore, Cs ~ 1 IF and a few micron long mechanical resonator 
gives a coupling ratio C m /Cs ~ 0.1. 

With all the capacitance and inductance values in hand, we are now ready to 
compare the coupling strengths gc and gi- We have: 



M m t / 2($o/2tt) y(LJ) 
9l = u w — \I — ~ 0.01u w (52) 



and 



gc = ^ f 2e)2 '^ CJ) ~ °- 01 ^> (53) 

where note that the CPW is twice as long (21) as in the voltage node case (/) so that 
the 2nd mode and 1st mode frequencies are the same, allowing a direct comparison. 
Interestingly, although the inductance ratio M m /L is about an order of magnitude 
smaller than the capacitance ratio C m /Cs, the flux quantum addition energy term is 
about an order of magnitude larger than the charging energy term, bringing the two 
coupling strength terms ( |52i) and ( 1531) into line. The reason for the difference in the flux 
quantum addition energy and charging energy magnitudes can be seen more clearly by 
taking the ratio of the square root terms: 



'2 (<& /27r) 2 /(L w Q h/e" Rk_ ^ 

(2ey/(2CJ) An^LjC: 4ttZ ~ ' 1 > 

Note that this is essentially just the ratio of the universal quantum of resistance (i.e., 
von Klitzing constant Rk = 25.8 kO) to the CPW impedance and so with typical 
microwave device impedances in the several tens of Ohms together with the 4n factor 
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in the denominator, we expect this ratio to be typically much larger than one. Thus, 
with regards to the microwave mode-qubit coupling strengths, there is no advantage to 
be gained by implementing one qubit scheme over the other. 

Comparing the qubit-mechanical oscillator coupling strengths (I34I) and flSJ) , we 

have 

Ac , VC (C m /Cy) (Vdc/Idc) , vin - 3 (W4) 

— = An — — — — ~ 5 x 10 — — - — , (55) 

Al Vl (M m /L) R K 1 fi ' v ; 

where we have assumed geometrical factors VcVl ~ 1- Thus, applying a CPW 
current bias Id c ~ 5 mA in the flux qubit device gives the same order mechanical 
coupling strength as applying a CPW voltage bias Vdc ~ 1 V in the charge qubit 
device. For a 5 micron long mechanical resonator with fundamental flexural frequency 
u) m j (27r) = 50 MHz, zeropoint displacement x zp = 10~ 14 m, and mechanical beam-centre 
conductor gap d = 100 nm, we have for the absolute, charge qubit-mechanical coupling 
strength 

Ac %zp C m eVdc „ 1 Vdc 
Kc = — Vc ~ U.l — . 

fojJm d Cy, fruj m IV 

Therefore, a voltage bias V^ c ~ 1 V provides sufficiently strong charge qubit-mechanical 
oscillator coupling strength kq< for the scheme of reference [I]. 

The V ac (t) and I ac (t) terms in Hamiltonians (I35p and (j46p are used to drive the 
mechanical oscillator into a non-zero displacement amplitude state, as required by the 
scheme of the companion paper pQ to rapidly generate a mechanical superposition state. 
Working with Hamiltonian ( 1351) and assuming a sinusoidal voltage drive that is resonant 
with the mechanical frequency, the steady state mechanical amplitude x Q m in units of 
the zeropoint displacement is 

t° t n r it/, it/ 

l"0| = = VC—T-^ Z Qmi (56) 

x zp d C s hw m 

where V® c is the ac voltage drive amplitude and Q m is the mechanical quality factor. 
For the same mechanical oscillator parameters as above and quality factor Q m = 10 4 
and dc voltage bias Vd c = 1 V, we have 

Kl ~ 10 7 ^|. (57) 

Thus, according to the mechanical superposition generation scheme P, ac voltage 
drive amplitudes of no more than a few microvolts are sufficient to obtain the required 
displacement amplitudes. Such ac voltage amplitudes do not significantly perturb the 
dc bias conditions on the CPB, since they are considerably smaller than the voltage bias 
difference 2e/C m ~ 0.5 mV separating neighbouring CPB charge degeneracy points. 

Alternatively, working with the Hamiltonian (1461) and assuming a sinusoidal current 
drive that is resonant with the mechanical frequency, the steady state mechanical 
amplitude in units of the zeropoint displacement is 

r° r M M \L \I° 
|«o| - — - Vl—, — z z Vm, (5SJ 
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where J° c is the ac current drive amplitude. Using the same mechanical resonator 
parameters as above and a dc current bias I^ c = 5 mA (for which the flux 
qubit-mechanical oscillator coupling strength Al coincides with charge qubit-oscillator 
strength Ac - see above), we have 

Kl ~ 10 6 ^|. (59) 

Thus, ac current drive amplitudes of no more than a few tens of microamps are sufficient 
to obtain the required displacement amplitudes for the mechanical superposition 
generation scheme [I]. Such ac current amplitudes do not significantly perturb the 
dc bias conditions on the flux qubit, since they are considerably smaller than the 
current bias difference 2<3>o/M m « 1 mA separating neighbouring symmetric double 
well configurations. 

The hxg(a w + a+)(a m + a+)cr x interaction terms in Hamiltonians fl35l) and f|46l) can 
in principle be employed to passively cool the mechanical oscillator (i.e., flexural mode) 
by driving the microwave mode off-resonance on the red-detuned side |10j . Provided 
g/\ui a — uj w \ <C 1, the described interaction term can be approximately replaced with 
the dispersive interaction: 

hxg(a w + a+)(a m + a^)a x — > - ;a+a w (a m + a+)cx 2 . (60) 

If the qubit is prepared and subsequently remains in its ground state during the driving 
interval, then we obtain the familiar ponderomotive interaction between the microwave 
and mechanical oscillator modes. For the above considered parameter values, we have 
for the dispersive coupling strength in either qubit case: 

~ Xq2 « ^ ~ 10- 9 . (61) 



Interestingly, this upper bound on the achievable coupling strength is still orders of 
magnitude smaller than what is possible in a related scheme where the SQUID with 
mechanically compliant loop segment is instead imbedded within the planar waveguide, 
interrupting the centre stripline conductor [TO]. Applying the theoretical analysis of 
reference [10] assuming a dispersive coupling strength ~ 10~ 10 , microwave mode of 
5 GHz and quality factor Q w = 10 4 , and mechanical flexural mode of 50 MHz (thus 
we are in the good cavity limit: uj m /j w 3> 1 [22]) and quality factor Q m = 10 4 , we 
find that pump drive powers of magnitude ~ 10 mW are required in order to cool 
the flexural mode down to a few energy quanta starting from a temperature of, say, 
100 mK. Such drive power strengths are unrealistic given the need to substantially 
attenuate the thermal Johnson noise in the original GHz drive signal entering the 
RF line at room temperature [H]. Thus, the dispersive coupling in the presently 
considered scheme is not strong enough in order to achieve significant passive cooling. 
Nevertheless, as emphasized in the companion paper pQ, there is no upper thermal 
threshold to verifying mechanical superposition states, although for the considered 
device parameters and temperatures higher than a few tens of mK, we expect it to 
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become difficult to tease out the recoherences in an echo experiment that are the 
signature of mechanical superposition states. As long as the mechanical resonator 
temperature can be maintained to within a few tens of mK (i.e., dilution fridge 
temperatures) during the drive and subsequent qubit control stages, there should be 
no problem observing mechanical superposition states; cooling is not essential for the 
scheme. 

5. Conclusion 

We have described an implementation of the nanomechanical quantum superposition 
generation and verification scheme presented in the companion reference pp. The 
implementation is based on the circuit QED set-up with the incorporation of a 
mechanical degree of freedom formed out of a suspended SQUID loop segment located 
opposite the centre conductor of a coplanar waveguide. Two qubit realizations were 
investigated, namely the capacitively coupled charge qubit with voltage-biased qubit- 
mechanical resonator coupling, and the inductively coupled flux qubit with current- 
biased qubit-mechanical resonator coupling. Both qubit realizations were found to have 
comparable and feasible resonator coupling strengths that are adequate for the quantum 
superposition generation scheme pQ. Ultimately, less predictable qubit properties such 
as their coherence times will decide which implementation is the more promising of the 
two. Nevertheless, the present considerations suggest that both implementations are 
equally worth pursuing in experiment. 
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Appendix A. Derivation of the single microwave mode-SQUID equations 

The most general solution to the wave equation (Tj[|) for <p(z, t) that satisfies the boundary 
conditions I w (±l/2,t) = takes the following Fourier integral form for — 1/2 < z < 0: 



with (p + (z,t) = ±0 (— z,t), < z < 1/2, for a voltage antinode (node) solution. 
Substituting Eq. (TO) evaluated at z = into Eqs. (TL2l). (jig). (E21). and (j23]l can serve 
as a starting point for analyzing the classical dynamics of the coupled CPW-SQUID 




dkA(k) cos(kvt) [cos(fcz + 0(h)) + cos(kz + kl - 6{k))\ 




(A.l) 
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system. Because of the nonlinear "force" terms in the SQUID part of the equations of 
motion ([TBI and fl22|) for 7± , the coupled dynamics is expected to be quite nontrivial 
and rich in general. On the other hand, if the force terms had been linear, then a 
simpler normal mode analysis could be used with only two distinct k values in (1A.1I) 
necessary for constructing the normal mode solutions. Furthermore, provided the two 
normal mode k values differ by much less than Z _1 [so that the dispersive differences in 
the ^-dependent terms of flA.ll) can be neglected] , then the wave equation describing 
the CPW can be approximately replaced by a much simpler, single mode harmonic 
oscillator equation with linear coupling to the linearized SQUID. In other words, the 
CPW-linearized SQUID system can be modelled simply as two coupled oscillators. 
Nevertheless, under conditions of weak coupling between the CPW and SQUID, i.e., 
C m <C C-z,C w l and M m <C L,L w l, it may still be possible to approximately replace 
the CPW with a single microwave frequency harmonic oscillator bilinearly coupled to 
the nonlinear SQUID. We shall assume this is the case and will now proceed in this 
appendix to derive the approximate single harmonic oscillator-SQUID coupled equations 
of motion for the voltage antinode, capacitively coupled case. We stress, however, that 
there still remains the issue of establishing the full domain of accuracy of the solutions 
to these approximate equations by comparing with the classical solutions assuming the 
full multimode Fourier integral form (1A.1I) . 

Restricting to single mode forms of equation ( 1A.1I) . we have 

0~O, t) = 2A cos[k(z + 1/2)] cos(kvt), (A.2) 

where we have fixed the phase in equation ( 1A.1I) to be 8(k) = kl/2. Defining 
0o = 2A cos(kl/2), the phase amplitude at the SQUID location z = 0, we have 

0- (0, t) = 0o cos(kvt) = (p{t) (A.3) 

and 

,90- OM) 



dz 



= -ktan(kl/2)(f) cos(kvt) = -ktm(kl/2)<p(t), (A.4) 

=o 

Substituting equations ( 1A.3I) and ( 1A.4j) into equation ( fl2i) and linearized equation ( {TBI) 



and expressing the latter two using more concise notation, we have 

7 + + «7 + = e0 (A. 5) 

and 

i0-/^tan(H/2)0=i 7+ , (A.6) 

where = -w 2 = -(kv) 2 (f), a = 4vr/ c cos(7r$ ext /$ )/(*oCs), P = [C^L^ 1 and 
e = C m /Cj] ^ 1. From the form of Eqs. ( 1A.5I) and ( 1A.6I) . we make a guess that the 
coupled oscillator equations of motion that yield approximately the same solutions take 
the form: 

7+ + a 7+ = 6e0 

+ c = de 7+ . (A.7) 
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However, because of the tan(&£/2) term in Eq. (1A.6[) . it will only be possible to match 
the solutions to those of the coupled oscillator equations with to-be-determined constant 
coefficients a, b, c and d, provided that the two normal mode frequencies are sufficiently 
close to each other so that dispersive differences can be neglected. 

In order to determine the oscillator equation coefficients in terms of a, {3 and e, 
we match the normal mode frequencies of Eqs. (1A.5|) and (1A.6I) to those of Eq. (1A.7|) . 
Solving first for the normal frequencies of Eq. (I A. 71) to second order in e, we obtain 

c bd 



UJ+ = k + v = y/o~ ( 1 + l^T_ e ^ ) ( A . 

2 a 



u - = k~v = yfc ( 1 - ^J^—e 2 ) . (A.9) 



and 



2a 

The secular equation following from Eqs. (1A.5I) and (1A.6I) is 



10 



2 



(3uJ ( Ujl 



-e H tan 



■1 



+ ^ = 0, (A.10) 



2 v \^vj\ 2 

where we have used k = uo/v. Solving Eq. (lA.lOj) for the smallest mode u to second 
order in e, we obtain: 

u+ = ^+™cot(^)e 2 (A.li; 

and 



4/3 V 2v 



v-^fa _^_(2W0 J »VW 2 (A12) 



/V ^ / a-(27rv//) 2 
Comparing Eqs. (1A. 1 1[) and 0A.120 with Eqs. (jA.80 and (IA.90 . we can clearly see that 
the CPW-linearized SQUID mode frequency uj + has a different form than the uj + arising 
from the coupled oscillator equations (I A. 7|) . This is because the boundary conditions 
on the CPW result in mode solutions with nonlinear dispersion, as signified by the 
presence of trigonometric terms. However, if we assume that the modes u + and uo~ are 
sufficiently close to each other, then we can expand the cotangent term in the mode 
frequency difference. Eq. ( 1A.11I) then becomes approximately 

+ (2nv/l) 3 v*/(/3l) 2 

Validity of the second order in e expansion together with this expansion in the frequency 
difference demands the following condition: 

m r 1 < \k + - k~\ <Z _1 . (A.14) 



Cy,C w I 

The CPW-linearized SQUID mode frequencies now take the same form as those for the 
coupled harmonic oscillators. Matching to determine the oscillator equation coefficients, 
we have: a = a, c = (2tcv/1) 2 (1 — 2v 2 e/((3l)), 6=1, and d = 2v 2 /(j3l). Substituting 
these values into Eq. (1A.7j) and restoring the nonlinear 7 + force term, we finally 
obtain the desired coupled single oscillator-SQUID equations of motion ( IT41) and ( TL5l) . 
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From Eq. (1A.14[) . we assume these equations are a good approximation to the full 
equations provided the difference in the characteristic frequencies for the SQUID and 
CPW dynamics, Su, satisfies 

°™ V < \Su\ < j. (A.15) 



Cj]C w l I I 
The derivation of the single oscillator-SQUID equations for the voltage node, 
inductively coupled case run along similar lines, yielding equations and (123]) . We 
assume these equations are a good approximation to the full equations provided the 
difference in the characteristic frequencies for the SQUID and CPW dynamics, doo, 
satisfies 

CJ ( M m \ 2 v v 

< \olj\ < -. (A.16) 



Cj V L J I 1 1 / 
Appendix B. Derivation of the qubit Langevin equations 

In this appendix, we derive the qubit Langevin equations. This operator formulation 
provides a complementary approach for addressing the dissipative qubit dynamics to 
that of the more commonly employed Fokker-Planck, master equation approach which 
solves for the time-dependence of the qubit density matrix pQ. The Langevin equation 
formulation of the qubit dynamics naturally belongs to the 'in-out' quantum optics 
approach to quantum measurement and control [23] and can be straightforwardly applied 
to analyzing circuit QED set-ups [10] . The qubit baths are modeled as a dense spectrum 
of oscillator modes and a rotating wave approximation is made for the assumed weak 
qubit-bath couplings: 



# q b-bath/ft= J dujuja{{uj)ai{uj) + / duuaT, (u)a 2 {uj) 
+ J duj[K*(uj)af(uj)a~ + K 1 (u)a + ai(u)) 

+ J dw[KZ(u)c£(u)<T, + K 2 (u})a z a2(w)], (B.l) 

where a + and a~ are the qubit (spin) raising and lowering operators [o x = a + + a~) and 
the ai, a 2 are independent oscillator bath modes, the first coupling K\ induces qubit 
decay through spin flips and where the second coupling K 2 causes pure dephasing. The 
equations for the bath modes are 

dai . -is* - 
dt 1 

= — %uja\ — iK 2 a z . (B.2) 

dt 

For the qubit, it suffices to consider the equations for a z and a + : 

-p. = - % \ 0z ,H\ -2i du{K 1 a + a 1 - K^at<r~) (B.3) 



and 



da + i 



dt 



^[a+, H] - i J duK; a fa z + 2i j duj(K*a+a + + K 2 a + a 2 ), (B.4) 
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where H is the charge qubit Hamiltonian (1351) or flux qubit Hamiltonian (j46p . 
Integrating the bath equations of motion (1B.2|) . we have 



ai (u,t)= -iK*(u) / dt'e- luj{t - t ' ) o--(t') + e- iuj{t - to) a 1 (LU,t ) 

J to 

a 2 (tu,t) = -iK*{u) [ dt'e-^-^a^t') + e- luj{t - to) a 2 (uj,t ). (B.5) 



to 



Substituting these solutions into equations ( IB. 31) and ( IB. 41) for a z and <r + , we have 



da z 
dt 



h 



dt' 



to 



a 



(t)a-(t') / dulK^e-^-^ 



+a + (t')o--(t) J du\K 1 \ 2 e +i ^ t -^ 
-2ia + {t) J duK^-^-^a^uo.to) 
+ 2i [ dLuKle +iu > {t - tQ) at{uo,tv)o-{t) 



(B.6) 



and 



da + 
~dT 



l -\a + ,H}+ f dt'a + {t')a z {t) [ dulK^e^-^ 
n Jt J 

i I duKle iuj(t - to) af(u,t )a z (t) 



-2 dt' 



'to 



-a 



a z (t')a + (t) J diulK^e^-^ 

~{t)a z {t') J d^K^e- 1 ^-^ 
+ 2i J duoK* 2 e iw{t - ta) a+(uj^)o + {t) 
+ 2ia + {t) J duK 2 e- iuj(t - to) a 2 (tj,t ). 

We now make the so-called first Markov approximation 
dependences of the bath couplings Ki to obtain 



(B.7) 

neglecting the frequency 



^ = -~[a z , H] - 7i (a z + 1) - 2zV7i"(e^V + < - e~'^a[ n+ a-) (B.8) 



and 



~dt 



l -[o-+,H]-( ll /2 + lv )o- 



- i^lie-'^a^(7 z + i^/2^(e- ilp a 2 n+ a + + e^a + a 2 n ), (B.9) 

where we have reparametrized Ki = \pfyj '(27r)e l<jil in terms of the decay rate 71, phase 
0x and K 2 = ^Ay^/ {^)e 19 in terms of the pure dephasing rate 7^ and phase <p. We 
have also used the identities a + o~~ = \(o~ z + /) and cr z a + = —a + a z = a + . The 'in' bath 
operators are defined as a- n (t) = (2tt)^ 1 ^ 2 J dcue~ luJ ^~ toS, ai(uj,t ). 
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